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ESTIMATION  OF  THE  SEA  SURFACE  RADAR 
CROSS  SECTION  AT  HF  FROM  SECOND-ORDER 
DOPPLER  SPECTRUM  CHARACTERISTICS 


INTRODUCTION 

The  radar  scattering  coefficient  of  the  sea  surface  at  HF,  <r°,  is  generally  difficult  to  obtain  by  a 
straightforward  well-calibrated  radar  measurement  because  several  of  the  variables  used  in  the  radar 
equation  are  not  well  known.  For  either  surface  wave  or  skywave  illumination  of  the  sea  surface,  the 
gain  to  be  used  in  the  radar  equation  is  difficult  to  estimate  because  it  can  be  a  complex  function  of 
local  ground,  soil  moisture,  ground  screen  quality,  coupling  to  the  sea  surface,  tidal  effects,  and  so  on. 
For  the  case  of  a  skywave  measurement  the  problem  is  further  compounded  by  the  variation  of  the  gain 
with  elevation  angle,  as  well  as  the  ionospheric  propagation  losses  the  signal  suffers  within  the  iono¬ 
sphere  in  its  two  refractive  passes. 

Knowledge  of  <r°  of  the  sea  by  a  means  independent  of  the  radar  equation  would  allow  one  to 
determine  the  unknowns  discussed  above,  as  well  as  prove  useful  for  comparison  with  targets  which 
might  be  detected  in  the  Doppler  spectrum  one  is  measuring.  This  would  allow  one,  for  example,  to 
attempt  to  apply  target  identification  techniques  to  the  target  of  interest.  These  have  been  developed 
over  the  past  several  years  by  Ksiensky  and  others  at  the  Ohio  State  University  Electroscience  Labora¬ 
tory  and  are  particularly  applicable  to  the  HF  band  of  radar  frequencies  [1],  These  techniques  depend 
upon  a  multi-parameter  measurement  of  a  target  and  a  comparison  of  that  measurement  with  a  catalog 
of  values  of  cross  section  as  a  function  of  the  parameters  available.  The  parameters  in  question  typi¬ 
cally  include  radar  frequency,  polarization  and  phase.  For  ionospheric  propagation,  radar  frequency  can 
be  varied  over  a  reasonable  bandwidth  for  coverage  of  a  given  illuminated  area,  of  the  order  of  6  MHz. 
For  the  elliptically  polarized  wave  which  exits  the  ionosphere,  a  measurement  of  the  difference  between 
the  minima  and  maxima  in  signal  amplitude  can  be  used  to  draw  some  conclusions  about  the  ratio  of 
horizontal  to  vertical  radar  cross  section  of  a  target.  Phase  information  in  an  absolute  or  even  relative 
sense  between  two  radar  frequencies  is  not  known  for  ionospheric  propagation.  Hence,  with  some 
information  about  two  of  the  three  parameters  available,  target  identification  with  skywave  HF  radar 
might  be  possible  with  a  satisfactory  calibration  cross  section  for  comparison,  such  as  that  of  the  sea 
surface. 

Another  application  of  such  cross  section  information  is  that  of  remote  sensing  of  the  sea  surface. 
It  is  well  known  that,  to  first  order,  the  cross  section  of  the  sea  surface  is  proportional  to  the  com¬ 
ponents  of  the  ocean  wave  directional  energy  spectrum  which  are  travelling  toward  and  away  from  the 
radar  bearing.  Hence,  assuming  one  knows  or  can  estimate  the  directional  spreading  of  the  wave  spec¬ 
trum,  one  can  then  determine  how  highly  the  sea  surface  is  developed.  Coupled  with  knowledge  of  the 
longest  wavelengths  excited,  and  an  estimate  of  wave  spreading  with  angle,  one  can  get  a  good  estimate 
of  RMS  waveheight. 

In  the  following  sections  we  test  a  simple  scheme  for  estimating  the  radar  cross  section  of  the  sea 
usmg  a  measurement  of  relative  amplitudes  of  the  approach  and  recede  Bragg  lines,  in  addition  to  the 
relative  amplitude  of  second-order  energy  around  zero  Doppler  frequency.  By  using  a  nomograph 
based  on  a  theoretical  model  for  the  radar  Doppler  spectrum  and  these  measurements,  a  value  of  «r°  is 
obtained.  Some  comparisons  are  made  with  surface  wave  Doppler  spectra  for  which  surface  truth  ocean 
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wave  spectra  are  available.  The  assumptions  made  in  the  theory  are  discussed  in  light  of  the  compari¬ 
son.  The  technique  is  particularly  useful  in  the  sense  that  it  can  deal  with  cases  in  which  the  actual  pre¬ 
vailing  ocean  wave  spectrum  may  have  severe  departures  from  idealized  spectra  (such  as  the  Pierson- 
Moskowitz  or  JONSWOP  models),  which  is  in  fact  the  case  for  several  of  the  data  comparisons.  The 
key  to  this  technique  is  that  both  radar  measurements  used  are  a  function  of  wave  components  in 
nearly  the  same  region  of  frequencies  of  the  wave  spectrum,  and  are  not  dependent  upon  the  long 
wavelength  regions  which  contribute  heavily  to  the  total  energy  of  the  ocean  wave  spectrum,  and  in 
which  much  of  the  amplitude  variability  may  occur.  It  is  noted  that  the  empirical  technique  developed 
in  Ahearn  et  al.  [2]  is  similar  to  the  one  developed  here,  except  for  directional  effects  and  an  absolute 
calibration. 

FIRST-ORDER  BRAGG  SCATTER 

The  integrated  or  total  radar  cross  section  measured  by  a  Doppler  radar  is  defined  by  the  equation: 

cr°  =  1/2  f  <r°(<o)  dot  ”  f  <r°(ci»)  d<o,  (1) 

*/  —  OO  0 

where  w  is  the  radio  frequency  in  radians.  In  operating  a  pulse-Doppler  radar  with  a  given  receiver 
band  pass  filter,  the  integrand  of  Eq.  (I)  is  multiplied  by  the  filter  function  describing  the  filter,  which 
effectively  changes  the  limits  of  integration  to  the  upper  and  lower  frequencies  of  the  pass  band.  Since 
the  radar  is  operated  at  some  pulse  repetition  frequency,  or  PRF,  the  sampling  theorem  effectively  folds 
the  energy  contained  in  the  pass  band  filter  into  a  narrower  bane  of  frequency  equal  to  that  of  the  radar 
PRF,  generally  tens  of  hertz  for  an  HF  radar.  The  PRF  is  typically  chosen  so  as  to  contain  all  of  the 
scattered  radar  energy  of  interest  in  an  unambiguous  manner  in  Doppler.  The  result  is  that  the  limits 
of  integration  of  Eq.  (1)  can  be  further  reduced  to  that  of  plus  and  minus  one  half  of  the  PRF.  The 
variable  of  integration  is  now  Doppler  frequency,  <oD.  (Unless  further  filtering  is  done,  the  noise  power 
rises  within  the  PRF  bandwidth  by  a  factor  equal  to  the  ratio  of  receiver  filter  pass  bandwidth  divided 
by  the  PRF,  assuming  it  is  white  noise  across  the  pass  band.  We  are  not  concerned  with  noise  here  so 
we  shall  ignore  this  effect  hereafter.) 

The  term  within  the  kernel  of  Eq.  (1)  is  referred  to  as  the  scattering  cross  section  per  unit  fre¬ 
quency,  or  scattering  cross  section,  as  opposed  to  the  total  or  integrated  cross  section  on  the  left-hand 
side.  For  sea  scatter  at  HF,  it  can  be  expanded  into  first,  second,  and  higher  order  contributions  in  the 
following  manner: 

cr°  (wp)  —  o-ft)  (wp)  +  <r^2)  (&»p)  + -  (2) 

where  the  two  terms  on  the  right  side  are  the  first-  and  second-order  contributions  to  the  radar  cross 
section. 

The  relationship  of  the  first-order  approach-recede  Bragg  line  contributions  to  the  radar  Doppler 
spectrum  with  the  wind  driven  ocean  wave  spectrum  is: 

<rfi)(0,<iiD)  —  8  it  Kg  [S(Kb,  0)8(a>p  +  <oB)  +  S(KB,  9  +  n)  S(wD  —  o>B)\,  (3) 

where  S(KB,  9)  is  the  directional  ocean  wave  spectrum;  9  is  the  angle  between  the  transmitted  radar 
vector  and  the  wind  direction;  KB  is  the  Bragg  resonant  ocean  wave  number,  related  to  the  radar  wave 
number,  k,  by  KB  -  2k  sin  <t> ;  <t>  is  the  angle  between  the  normal  and  incident  (and  scattered)  radio 
wave;  and  <uB  —  fifl  —  y/gKB  is  the  Bragg  Doppler  frequency,  equal  to  the  ocean  wave  Bragg  fre¬ 
quency.  Barrick  [3]  derived  Eq.  (3)  for  HF  applications,  and  a  surface  wave  experiment  by  Teague  et 
al.  [4],  verified  the  theory  to  within  3  dB  accuracy.  That  experiment  made  use  of  measurement  of 
energy  backscattered  from  the  sea  made  simultaneously  with  energy  received  from  the  transmit  site,  a 
LORAN  C  station.  Hence,  antenna  gains  of  both  the  transmitter  and  receiver  cancel  one  another  when 
taking  a  ratio  of  radar  equations  for  both  paths  (received  power  for  the  direct  path  divided  by  received 
power  for  sea  backscatter) .  The  radar  data  agreed  to  within  3  dB  with  a  buoy  measurement  of  the 
omnidirectional  spectrum  at  the  Bragg  resonant  frequency,  0.14  Hz. 
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One  can  write  the  ocean  wave  directional  power  spectrum  as: 

S(K)  -  S(K.O)  -  SUO  F(9),  (4) 

where  the  approximation  has  been  made  that  the  angular  dependence  of  the  spreading  or  ocean  waves 
is  independent  of  wave  number.  The  first  term  is  the  so-called  omnidirectional  energy  spectrum,  since 
F(8)  integrated  over  all  angles  gives  unity.  A  typical  form  for  the  omnidirectional  spectrum  is  the 
Pierson-Moskowitz  spectrum,  [5]: 

S(n)<m  =  ag 2  H-5  exp  (-  big/  (l  w)4)  dCl  (5a) 

or,  in  terms  of  wave  number, 

SiK)  dK  -  (a/2)  K  *  exp  i-big/Kw2)2)  dK,  (5b) 

where  w  is  the  wind  speed  at  19.4  ft  above  the  sea  surface;  b  is  a  constant,  0.74;  fl  is  the  wave  fre¬ 
quency  in  radians;  and  a  is  typically  of  the  order  of  0.0081  for  long  fetches,  as  proposed  by  Phillips,  but 
which  has  been  shown  to  be  a  function  of  fetch.  (See,  for  example,  Hasslemann  et  al.,  [61.)  A  plot  of 
a  series  of  ocean  wave  spectra  is  shown  in  Pig.  1  as  a  function  of  wind  speed  generating  the  seas,  and 
also  for  two  values  of  the  variable,  a.  The  effect  of  changing  a  is  to  raise  or  lower  the  asymptotic  H~s 
line. 


Fig.  I  —  Two  families  of  curves  representing  ocean  wave  spectra  given  by  Eq. 
(Sb)  in  the  text  are  shown.  Each  family  is  parametric  in  the  wind  speed,  with 
the  low  frequency  amplitudes  growing  with  the  wind.  The  difference  between 
the  two  families  is  the  constant,  a,  the  upper  family  determined  using  a  value 
of  0.0081,  while  the  lower  family  is  determined  by  an  a  of  0.00081.  The  deci¬ 
bel  difference  between  the  two  asymptotes  at  the  high  frequency  end  is  seen  to 
be  ten  dB  as  a  result  of  the  factor  of  ten  decrease  in  a. 


Since  the  radar  cross  section  is  determined  essentially  by  the  sum  of  the  first-order  contributions 
to  the  Doppler  spectrum,  and  these  terms  are  proportional  to  the  amplitude  of  the  ocean  wave  spec¬ 
trum  at  the  Bragg  wavelength,  the  radar  cross  section  will  vary  linearly  with  the  parameter  a: 

< r°(0,toD )  =  4 it  a  (F(9)bi<i>D  +  ®r)  +  F{9  +  rr)6(a»o  —  &>g)).  (6) 

In  the  above  we  have  set  the  exponential  factor  of  Eq.  (5b)  to  one,  since  it  affects  the  spectrum  near 
the  low  frequency  falloff,  but  not  near  KB  (5  to  15-m  ocean  wavelengths). 
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When  operating  in  a  normal  monostatic  mode,  if  one  steers  a  very  narrow  beam  radar  operating 
line  of  sight  or  in  the  surface  wave  mode,  one  will  observe  the  maximum  Bragg  line  of  the  approach- 
recede  pair  to  reach  its  maximum  amplitude  in  the  direction  that  is  parallel  to  the  wind.  The  weaker  of 
the  two  will  likewise  reach  its  minimum  value  at  this  bearing.  The  ratio  of  the  two,  p,  will  thereby 
reach  its  maximum  value.  At  right  angles  to  the  wind  direction,  this  decibel  ratio  goes  to  zero,  as  the 
maximum  Bragg  peak  has  decreased  and  the  minimum  has  increased  such  that  the  two  are  equal.  This 
behavior  with  angle  reflects  the  fact  that  the  variations  in  p  are  strictly  due  to  the  variation  of  the 
spreading  function  of  Eq.  (4)  with  angle.  Since  the  omnidirectional  spectrum  is  the  same  over  the  area 
covered  by  the  radar  (very  short  ranges  are  assumed),  the  only  variation  in  the  left  hand  side  of  Eq.  (4) 
is  due  to  F(e).  Given  a  model  for  the  spreading  function,  one  can  use  this  technique  to  estimate  the 
angle  that  the  radar  bearing  makes  with  the  wind  driven  waves  for  a  single  radar  bearing.  This  tech¬ 
nique  has  been  implemented  in  the  skywave  mode  of  operation  to  derive  wind  direction  fields  over 
large  expanses  of  the  ocean  [7]  and  to  locate  hurricane  centers  and  define  fronts  using  very  narrow 
beam  antennas  [8]. 

Although  this  technique  allows  the  estimation  of  direction,  the  problem  of  estimating  the  first 
term  of  Eq.  (4)  from  first-order  Bragg  line  measurements  used  with  the  radar  equation  alone  remains 
difficult  because  of  the  antenna  gain  uncertainties.  Additional  information  is  available  in  the  second- 
order  contribution,  however,  which  can  allow  one  to  estimate  the  amplitude  of  the  parameter,  a,  and, 
hence,  allow  a  similar  estimation  of  a0.  First,  an  introduction  into  the  second-order  contribution  to  the 
Doppler  spectrum  is  necessary. 

SECOND-ORDER  BRAGG  SCATTER 

In  the  following,  we  shall  give  an  intuitive  development  of  the  simple  corner  reflector  model  used 
to  describe  second-order  scatter.  This  is  the  basis  of  the  work  by  Trizna  et  al.  [9)  and  differs  from  work 
by  Barrick  [10],  who  has  considered  both  electromagnetic  and  hydrodynamic  contributions  to  the 
second-order  structure  of  the  radar  Doppler  spectrum.  Valenzuela  [11]  has  considered  both  of  these 
contributions  as  well,  but  his  results  are  more  general,  since  they  were  derived  for  microwave  frequen¬ 
cies,  and  are  valid  at  the  long  wavelength  limit  appropriate  to  HF  frequencies.  The  hydrodynamic 
resonant  term  of  Valenzuela  differs  from  that  of  Barrick  in  that  it  contains  a  viscous  damping  term,  v , 
which  is  very  important  for  the  microwave  radar  frequencies  examined  in  detail  by  Valenzuela,  and 
possibly  in  the  HF  regime  as  well.  This  point  is  an  area  of  current  investigation.  A  rigorous  develop¬ 
ment  of  the  corner  reflector  model  is  in  preparation  for  publication,  and  begins  with  the  integral  equa¬ 
tion  for  the  second-order  contribution  to  the  scattering  cross  section  developed  by  Valenzuela  (Trizna, 
in  preparation). 

Basic  Physical  Assumptions  for  Second-Order  Scatter 

For  the  second-order  electromagnetic  contribution,  which  is  the  major  term  near  zero  Doppler  for 
all  cases,  we  use  the  corner  reflector  model  as  outlined  in  Trizna  et  al.  [91.  The  electromagnetic  wave 
scatters  twice  from  all  possible  pair  combinations  of  ocean  waves  that  are  travelling  perpendicular  to 
one  another.  Each  individual  pair  produces  a  response  at  a  different  point  in  the  Doppler  spectrum 
which  is  determined  by  the  Doppler  shifts  produced  by  each  wave  travelling  at  its  appropriate  phase 
velocity.  The  electromagnetic  resonance  term  of  Barrick  allows  all  possible  combinations  of  ocean  wave 
vectors  to  contribute  in  double  scatter,  without  the  restriction  that  the  wave  pairs  be  perpendicular  to 
one  another.  Instead,  he  considers  the  first  wave  as  arbitrary  to  begin  with,  as  we  do,  then  considers  all 
possible  ocean  waves  for  the  second  scattered  radar  wave,  using  a  factor,  T,  which  is  resonant  for  the 
case  when  the  second  wave  is  perpendicular  to  the  first.  In  effect,  the  difference  between  the  two  tech¬ 
niques  is  that  he  considers  P  as  a  resonance  with  some  shape,  whereas  we,  in  effect,  evaluate  it  as  a  8- 
function  when  performing  the  integration. 
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The  problem  of  radio  scatter  from  physical  ocean  waves  is  considered  in  the  same  sense  as  one 
would  consider  radio  scatter  from  any  physical  object,  with  regard  to  the  calcula'ion  of  the  phase  shifts 
produced  by  the  nonzero  velocity  of  the  object.  Two  laws  of  physics  can  be  used  as  starting  points  for 
setting  up  the  calculation,  namely,  the  law  of  conservation  of  momentum,  and  the  law  of  conservation 
of  energy.  It  is  to  be  noted  that  these  are  in  fact  laws  of  physics  which  have  been  verified  and  used  in 
classical  mechanics  and  electromagnetism,  as  well  as  quantum  physics  and  special  and  general  relativity. 
They  are  starting  points,  for  example,  for  the  calculation  of  the  relativistic  interaction  of  photons  with 
matter,  and  degenerate  to  our  problem  in  the  classical  limit  (See,  for  example,  Leighton  [12],  Ch.  1). 

Double  Scatter  Wave  Numbers— Conservation  of  Momentum 

If  one  first  invokes  the  law  of  conservation  of  momentum  for  the  scatter,  then  the  following 
equation  results  for  the  total  momentum  conserved  for  the  double  scatter  of  the  electromagnetic  wave 
from  two  ocean  waves  with  wavenumbers,  K\  and  K 2: 

h  K\  +  h  K7  =  2  h  k,  (7) 

which  is  a  vector  equation,  where  * h '  is  Planck's  constant.  The  scatter  can  be  divided  into  two  indepen¬ 
dent  scatters,  with  an  intermediate  electromagnetic  wave,  with  wavenumber,  k'.  Momentum  must  be 
conserved  in  this  intermediate  scatter  as  well  ((13),  p.  145).  The  resulting  equation  describing  the  first 
scatter  is  given  by: 

hk  =  h  K]  +  hk'.  (8a) 

If  one  solves  the  above  equation  in  component  form,  one  sees  that  the  intermediate  scatter  must  reflect 
from  the  first  ocean  wave  at  an  angle  to  the  first  wave  number  vector  equal  and  opposite  to  the  incident 
angle.  One  can  now  consider  the  second  scatter,  with  the  incident  wave  now  the  intermediate  wave,  k'. 
A  similar  equation  results,  namely: 

hk'~  hK2  +  h  k".  (8b) 

For  the  radar  energy  to  scatter  back  to  the  radar,  k  "must  be  in  the  opposite  direction  as  k,  so 

k"=-k,  (9) 

where  the  '  denotes  the  unit  vector  defining  the  direction  for  each  wave  vector.  Solving  Eq.  (8b)  in 
component  form  again  requires  the  scattered  angle  be  equal  to  the  incident  angle.  The  choice  of  the 
direction  defined  by  Eq.  (9),  when  combined  with  (8a)  and  (8b),  sets  the  direction  of  K2  as  perpendic¬ 
ular  to  K |,  as  is  shown  in  Fig.  2.  That  is,  because  conservation  of  momentum  demands  that  radar 
waves  scatter  specularly  from  ocean  waves  in  the  plane  of  the  ocean  surface,  choosing  the  final  scatter 
direction  to  be  back  to  the  radar  requires  the  two  ocean  waves  participating  in  the  scatter  to  be  perpen¬ 
dicular  to  one  another.  Instead  of  this  intuitive  development  using  the  intermediate  scatter  concept, 
Eq.  (7)  can  be  used  alone,  in  conjunction  with  the  integral  equation  for  <r<°2)  (a»s,0),  as  a  delta  func¬ 
tion,  to  derive  the  closed  form  for  cr <02)  (u)p.9)  for  a  perfectly  conducting  sea  surface.  This  mathemati¬ 
cal  development  does  not  assume  that  the  intermediate  scatter  be  specular  in  the  plane  of  the  sea  sur¬ 
face,  but  it  follows  as  a  result  of  the  development. 

Doppler  Shift  Calculation— Conservation  of  Energy 

For  a  given  pair  of  ocean  waves  satisfying  Eq.  (7),  the  amount  of  Doppler  shift  from  the  transmit¬ 
ted  frequency  is  determined  by  the  conservation  of  energy  equation: 

(10) 

where  ft  -  h/ 2ir  and  where  wp  is  the  Doppler  shift  of  the  received  radar  energy  scattered  via  ocean 
waves  I  and  2.  Note  that  there  are  more  than  one  pair  of  ocean  wave  frequencies  which  can  satisfy  Eq. 
(10)  and  all  such  pairs  must  be  considered. 
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A:  *,  =  *2  =  "I 

K,  K, 


Fig.  2  —  The  double  scatter  mechanism  responsible  for  second  order 
Doppler  structure  is  shown  by  a  plan  view  as  seen  from  above  the 
surface.  Four  possible  combinations  are  indicated  for  one  value  of 
the  angle  between  K ,  and  *.  B  -  45  deg,  with  cases  B  and  C  giving 
a  zero  Doppler  contribution.  Cases  A  and  D  will  contribute  at 
±2 }/*fg,  respectively.  A  case  for  a  second  value  of  B  is  shown  in 
dashed  lines,  with  different  wavelengths  for  K\  and  K2.  yielding  a 
different  Doppler  shift. 


Determination  of  the  Scattered  Energy  per  Doppler  Filter 

The  calculation  of  expected  energy  per  Doppler  fijter  is  carried  out  in  a  systematic  manner,  run¬ 
ning  through  the  combination  of  wave  vectors,  K j  and  K2,  which  satisfy  Eq.  (7).  One  begins  with  the 
magnitude  of  K\  very  nearly  zero,  and  K2  very  nearly  equal  to  the  Bragg  wave  number.  Kg.  The 
angles  that  the  two  wave  vectors  make  relative  to  the  incident  electromagnetic  wave  vector  are  required 
to  satisfy  the  vector  equation  of  Eq,  (7),  in  component  form  along  arbitrarily  chosen  orthogonal  axes, 
say  the  y-axis  along  the  incident  electromagnetic  wave  vector  in  Fig.  2.  The  corresponding  ocean  wave 
frequencies  are  determined  using  the  dispersion  relation  for  linear  gravity  waves  relating  KB  and  fla, 
given  just  after  Eq.  (3).  The  Doppler  frequency  at  which  this  contribution  appears  in  the  Doppler  spec¬ 
trum  is  then  given  by  Eq.  (10).  At  each  Doppler  spectral  point,  for  a  given  Doppler  filter  width,  there 
is  a  corresponding  spread  in  angle  and  wave  number  dK}  and  dK2,  which  contribute  to  the  Doppler 
filter  considered.  For  a  specified  look  angle  of  the  radar  relative  to  the  wind  direction,  the  amplitude 
contribution  at  each  Doppler  spectral  point  is  then  determined  by  integrating  the  expression  for  the 
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second-order  scattering  coefficient,  o-fo,  over  Doppler  frequency  with  an  appropriate  Doppler  filter 
shape  tc  determine  the  scattered  energy  within  the  given  Doppler  filter: 

0)  —  26nk4  j"  dap  J*  dK i  |T p  S(K j)  5 (A" 2)  G  (u)q  —  <u^)  8 (<u£>  ±  f)  ]  ±  IT 2),  (11) 

where  T  is  given  in  [10,11],  and  G(aiD  -  rn'iy)  represents  the  spectral  response  of  the  FFT  filter,  typi¬ 
cally  a  sinc-function  if  no  time  weighting  is  done  on  the  FFT.  This  integral  is  further  considered  in  a 
future  publication  (Trizna,  in  preparation),  but  the  results  are  essentially  the  same  as  in  [9], 

Note  that  the  first-order  Bragg  line  energy  is  obtained  in  a  similar  manner,  integrating  Eq.  (6) 
over  Doppler  frequency  with  an  appropriate  Doppler  filter  shape.  Note  also  that  the  first-  and  second- 
order  contributions  behave  differently  with  different  coherent  integration  times,  or  Doppler  filter  line 
widths.  Because  the  first-order  Bragg  line  is  a  coherent  discrete  spectral  line  return  with  fixed  total 
energy,  the  response  will  increase  relative  to  the  noise  level  with  narrowing  Doppler  filter  width.  How¬ 
ever,  the  second-order  term  is  a  noise-like  continuum,  and  as  the  Doppler  filter  is  narrowed,  a  smaller 
spread  in  dKx  and  dK 2  is  included  for  a  given  Doppler  filter,  as  discussed  above.  This  excluded  energy 
appears  in  an  adjoining  Doppler  filter,  and  as  such,  behaves  as  a  noise  spectrum  rather  than  a  coherent 
narrow  spectral  return.  (As  the  coherent  integration  time  is  doubled  in  the  FFT  process  and  the  spec¬ 
tral  width  of  the  FFT  filter  is  halved,  the  number  of  spectral  filters  is  doubled  across  the  total  Doppler 
bandwidth.)  The  details  of  this  discussion  will  be  covered  in  the  future  publication  cited  above. 
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TECHNIQUE  OF  <r°  DETERMINATION 

The  determination  of  a-0  depends  upon  two  measurements,  that  of  the  ratio  of  approach/ recede 
Bragg  line  amplitudes,  p,  and  a  measure  of  the  amplitude  ratio  of  maximum  Bragg  line  to  zero  Doppler 
continuum,  £,  plus  a  consideration  of  the  coherent  integration  time  used  in  the  processing.  The  first 
measurement  allows  one  to  determine  the  angle  to  wind/wave  direction  the  radar  bearing  makes,  and 
thereby  allows  a  determination  of  F(9)  in  Eq.  (3).  The  second  measurement  allows  a  determination  of 
the  Phillips'  constant,  a,  of  Eq.  (5),  using  Eqs.  (6)  and  (11)  for  zero  Doppler  frequency.  With 
knowledge  of  a  and  9  the  radar  cross  section  is  thereby  determined  using  these  in  Eq.  (6).  A  final 
correction  is  also  made  for  the  integration  time,  or  Doppler  filter  width  used,  and  is  included  in  the 
final  results. 

Figure  3  shows  two  spectra,  illustrating  the  effect  of  a  different  value  of  a.  For  the  same  wind 
speed,  but  a  10  dB  decrease  in  a  (vis-a-vis  Fig.  1),  the  Doppler  spectrum  frequency  characteristics 
remain  quite  similar,  except  for  a  10  dB  decrease  in  the  first-order  Bragg  lines,  and  a  20  dB  decrease  in 
the  second-order  continuum  amplitude.  This  is  what  one  would  expect  for  once-  and  twice-scattered 
radar  waves,  each  scatter  being  proportional  to  a  resonant  ocean  wave  component,  vis-a-vis  Eqs.  (3) 
and  (11). 

Determination  of  the  Relevant  Ocean  Wave  Components 

The  relationship  between  measured  £  and  the  resulting  Phillips  constant,  a,  requires  knowledge  of 
which  components  are  being  sensed  by  the  radar  at  zero  Doppler.  From  Eq.  (10),  the  Doppler  fre¬ 
quency  is  zero  for  the  corner  reflector  case  of  the  magnitudes  of  the  two  ocean  wave  frequencies  being 
equal,  which  requires  that  the  magnitude  of  the  corresponding  wave  numbers  also  be  equal.  In  addi¬ 
tion,  if  one  writes  Eq.  (7)  in  component  form,  in  a  set  of  coordinates  discussed  earlier,  one  has: 

jj  K\  cos  0  +  s2  Kj  sin  ®  (12a) 

si  K]  sin  0  -  s2  K2  cos0  =  0  (12b) 
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DOPPLER  (f bragg*  DOPPLER  HbbaggI 

Fig.  3  —  Two  different  radar  Doppler  spectra  illustrate  the  effect  of  changing  the  constant  a  of  Eq  (4)  by  a  factor  of 
ten.  due  to  the  sea  surface  being  10  dB  down  from  traditional  saturation  values  in  the  high  frequency  regime.  The 
first-order  Bragg  lines  are  down  10  dB  and  the  second-order  continuum  is  down  20  dB.  ten  more  than  the  first  order, 
in  the  second  spectrum  compared  to  the  first.  These  might  be  the  result  of  wave  spectra  from  Fig.  I.  for  the  same 
wind  speed  and  the  two  ur's  differing  by  a  factor  of  ten 


where  0  is  the  angle  between  Kx  and  the  *-axis,  and  S]  and  s2  are  arbitrarily  chosen  as  +  or  -1.  The 
region  between  the  Bragg  lines  filled  when  sx  and  s2  are  chosen  with  opposite  signs,  and  to  the 
right/left  of  the  Bragg  lines  when  both  are  positive/negative.  The  only  0  which  satisfies  Eq.  (12b)  for 
Kx  =  K2  and  S|  =  -  s2  is  0  =  45°.  From  Eq.  (10),  for  zero  o>0,  fl,  —  12 2,  and  hence,  Kt  =  K2  for 
Ji  =  -  s2.  Substituting  this  value  into  Eq.  (12a)  then  gives  the  magnitude  of 

AT,  =  AT2  A:/(sin  45  °)  =  V2  A:.  (13) 

The  corresponding  ocean  wave  frequency  (using  the  gravity  wave  dispersion  relation,  fi  =  VgK,  and 
KB  =  2  k)  responsible  for  scattered  energy  as  zero  Doppler  frequency  is  found  to  be: 

n0=  n,  -  nfl/21/4  (H) 

-  0.841  ns. 

Hence,  the  ocean  wave  frequency  sensed  by  the  second-order  Doppler  energy  at  zero  Doppler  fre¬ 
quency  is  very  nearly  the  same  as  the  wave  frequency  sensed  by  the  first-order  Bragg  scatter,  and  is 
generally  far  from  the  peak  of  the  spectrum,  but  almost  certainly  in  the  region  where  the  wave  fre¬ 
quency  spectrum  behaves  as  fl~5.  The  directions  of  propagation  of  the  two  components,  such  that  the 
net  Doppler  frequency  is  zero,  is  determined  to  be  45°  to  the  right  and  left,  respectively,  of  the  radar 
pointing  direction,  with  one  ocean  wave  approaching  and  the  other  receding  from  the  radar.  (Note  that 
there  are  two  such  pairs,  cases  A  and  D  of  Fig.  2,  and  both  are  taken  account  of  in  the  radar  spectrum 
model.)  In  the  model  used  for  determination  of  (rfaffwp),  the  assumption  is  made  that  the  ocean  wave 
power  spectrum  does  vary  as  fl~5  near  the  Bragg  wave  frequency  in  amplitude,  with  the  angular  depen¬ 
dence  of  Eq.  (4)  varying  as: 

F(B)  -  —  (a  +  cos2  9),  for  0  <  |0|  <  w/2  (15) 

ir 

F(0)  —  2 a2/n  (a  +  cos2  0)_l,  for  ir/2  <  |0|  < 
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(See  [7]  for  a  further  description  of  this  model.)  Because  the  wave  number  determined  by  Eq.  (13)  is 
very  close  to  the  Bragg  wave  number,  the  error  in  the  Q -s  assumption  is  minimal.  A  greater  error  is 
probably  present  due  to  the  assumption  of  the  angular  spreading  used.  We  shall  discuss  this  in  a  later 
section. 

Since  the  Bragg  ocean  wave  frequency,  fls,  and  the  wave  frequency  contributing  to  the  zero 
Doppler  amplitude,  fi0,  are  nearly  equal  from  Eq.  (14),  both  are  unaffected  by  the  exponential  of  Eq. 
(5)  and  will  scale  in  a  similar  manner  with  changing  a.  That  is,  as  long  as  the  peak  of  the  wave  spec¬ 
trum  is  much  less  than  the  Bragg  ocean  wave  frequency,  slight  changes  in  the  position  of  the  peak  due 
to  different  wind  speeds  will  have  no  effect  on  the  wave  spectrum  amplitude  at  n  B  and  D0.  The  ampli¬ 
tude  at  these  frequencies  will  be  determined  only  by  the  parameter,  a,  (which  is  a  function  of  fetch) 
and  wind/ radar  angle,  so  long  as  the  wind  is  sufficiently  strong  to  drive  the  peak  of  the  spectrum 
far  from  the  Bragg  frequency,  so  that  the  exponential  term  of  Eq.  (5)  is  nearly  one. 

Results  in  Nomograph  Form 

The  model  for  the  Doppler  spectrum  by  [9]  was  used  to  generate  curves  of  £,  the  ratio  between 
the  strongest  Bragg  line  in  the  Doppler  spectrum,  and  the  zero-Doppler  continuum  amplitude,  for  vari¬ 
ous  values  of  the  parameter,  a.  The  results  are  shown  in  Fig.  4,  in  the  form  of  a  nomograph.  The 
horizontal  axis  is  wind  to  radar  angle.  The  positive  y-axis  is  p ,  the  ratio  of  approach  to  recede  Bragg 
line  power.  The  negative  y-axis  is  the  parameter,  £,  the  ratio  of  the  maximum  Bragg  line  power  to  the 
power  level  of  the  zero-Doppler  continuum,  which  is  a  function  of  wind-radar  angle.  The  curve  in  the 
lower  half  plane  labeled  with  A  equal  to  zero  is  the  locus  of  £  for  a  value  of  the  constant,  o,  equal  to 
0.0081,  historically  taken  as  the  Phillips  "saturation"  value  for  the  constant  a,  and  a  coherent  integra¬ 
tion  time  of  200  s. 


tig  4  —  A  nomograph  derived  from  theoretical  spectra  that  at-  ^ 
lows  one  to  calculate  the  number  of  dB.  A.  the  ocean  wave  " 
Bragg  wavenumber  (and  <r"l  is  down  from  the  Phillips  satura¬ 
tion  value.  00X1  (used  onlv  for  a  reference  valuer  With  the 
measured  value  of  ft.  the  approach-rccede  Bragg  line  ratio,  one 
determines  the  angle  to  wind  or  the  radar  bearing.  «.  using  the 
upper  half  plane  Measuring  C.  the  ratio  of  largest  Bragg  line  to 
zero  Doppler  second  order  amplitude,  one  draws  a  horizontal 
line  The  point  at  which  this  horizontal  line  intersects  the  verti¬ 
cal  line  drawn  down  from  the  estimated  radar-wind  angle,  «. 
determines  A.  the  number  of  dB  the  point  is  down  from  the 
theoretical  curve  for  the  specified  integration  time  This  value  is 
then  used  in  f.q  ’I  of  the  text  to  determine  <r"(«l 
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If  the  omnidirectional  sea  spectrum  is  down  6  dB  from  this  value,  then  the  first-order  Bragg  line 
amplitudes  will  also  be  down  6  dB  each  from  their  "saturated”  values.  The  ratio  of  the  Bragg  lines  will 
not  be  different,  however,  since  they  depend  only  upon  the  directional  factor  of  Eq.  (3).  Since  the 
zero-Doppler  amplitude  is  of  second  order  in  the  ocean  wave  spectrum,  it  will  be  down  12  dB  from  its 
value  for  the  "saturated"  sea.  The  ratio  £  will  decrease  by  the  difference  in  the  dB  change  of  the  first- 
order  Bragg  line  and  the  second-order  zero-Doppler  continuum,  i.e.,  6  dB.  A  second  curve  is  shown  in 
the  figure,  just  6  dB  below  the  first  with  the  same  variation  with  angle.  The  value  of  a  corresponding 
to  this  curve  is  therefore  0.0020  for  the  same  integration  time  as  the  first  curve.  This  second  curve  can 
also  be  considered  as  derived  for  a  equal  to  0.0081,  but  for  an  integration  time  four  times  less  than  for 
the  first  curve,  as  indicated  in  the  figure,  50  s.  (This  is  of  course  due  to  the  fact  that  the  Bragg  line  is  a 
coherent  signal,  while  the  continuum  is  broadband  and  hence  behaves  as  noise.  They  sum  like  T  and 
VT,  respectively,  where  T  is  the  coherent  integration  time.) 

To  use  the  curves,  one  takes  the  measurement  of  p,  draws  a  horizontal  line  until  it  intersects  the 
curve  in  the  upper  half  plane,  and  determines  a  radar-wind  angle  by  dropping  a  perpendicular  from  this 
intersection  to  the  x-axis.  One  continues  this  vertical  line  into  the  lower  half  plane.  One  then  takes  the 
value  of  £  measured  from  the  radar  spectrum  which  gave  p,  and  draws  a  horizontal  line  until  it  inter¬ 
sects  the  vertical  line  drawn  earlier.  The  number  of  dB  this  intersection  lies  below  the  theoretical  curve 
is  A,  the  number  of  dB  down  the  scattering  cross  section  of  the  sea  is  from  its  saturated  value  for  that 
radar  wind  angle.  This  value  must  yet  be  corrected  for  the  integration  period  one  used  in  processing 
the  radar  Doppler  spectrum  that  produced  p  and  £.  Since  the  curves  were  generated  for  a  200  s  integra¬ 
tion  time,  if  the  integration  time  used  is  half  this  value,  3  dB  should  be  added  to  A;  if  quadrupled,  6 
dB  is  added,  and  so  on.  In  the  case  of  ionospherically  propagated  signals,  particularly  via  F-layer 
modes,  one  should  either  choose  an  integration  time  such  that  the  Bragg  lines  just  fill  the  Doppler 
filter,  or  sum  up  the  energy  in  all  the  Doppler  filters  which  have  contributions  from  the  Bragg  line  for 
longer  integration  times.  The  simplest  course  is  to  use  50-s  integration  "times,  which  have  been  found 
optimum  for  very  good  ionospheric  conditions,  and  use  the  50-s  curve  in  the  figure.  If  the  ionospheric 
conditions  are  relatively  poor  with  a  50-s  value,  so  that  the  Bragg  line  contributes  to  several  adjacent 
Doppler  filters,  then  one  must  sum  all  of  the  energy  in  the  Bragg  lines,  possibly  encompassing  several 
Doppler  filters  in  calculating  £.  This  is  not  important  in  determining  p,  since  both  Bragg  lines  are 
affected  the  same  way,  but  must  be  done  in  measuring  £.  These  effects  are  discussed  in  greater  detail 
in  the  next  section. 

EFFECTS  OF  INTEGRATION  TIME,  RADAR  SIGNAL  BANDWIDTH, 

AND  EFFECTIVE  BANDWIDTH  DUE  TO  DEGRADATION 

OF  SEA  SURFACE  AND  IONOSPHERIC  COHERENCE 

There  are  several  limits  that  prevent  one  from  achieving  an  infinitely  narrow  Bragg  line  width  as 
one  increases  the  coherent  integration  time  indefinitely.  Barrick  indicated  in  his  work  that  the  Bragg 
line  has  zero  width  and  infinite  amplitude,  but  that  the  product  of  the  two  is  a  finite  value.  However, 
he  did  not  consider  the  fact  that  one  transmits  a  finite  bandwidth,  scatters  from  finite  patch  sizes  on  the 
sea  surface,  and  is  affected  by  coherence  lengths  within  the  ionosphere.  (His  calculation  assumed  an 
infinite  plane  wave  upon  a  perfectly  coherent  and  infinite  surface.)  We  consider  each  of  these  in 
increasing  order  of  importance  in  this  section. 

Doppler  processing  of  coherent  radar  data  is  equivalent  to  matched  filtering,  so  that  the  response 
from  the  filter  has  a  finite  width.  In  terms  of  the  FFT,  each  spectral  point  has  a  finite  width  equal  to 
the  total  FFT  bandwidth,  which  is  equal  to  the  pulse  repetition  frequency  (assuming  sum  and  quadra¬ 
ture  FFT  inputs),  divided  by  the  number  of  FIT  input  points.  If  the  backscatter  is  indeed  infinitely 
narrow  inherently,  the  measured  spectrum  will  equivalently  be  1/THz  wide,  where  T  is  the  coherent 
integration  time  used.  (Doubling  the  coherent  integration  time,  or  input  points  to  an  FFT.  is 
equivalent  to  halving  the  bandwidth  of  the  matched  filters  and  doubling  their  number.)  As  the  integra¬ 
tion  time  is  increased  indefinitely,  the  spectrum  would  not  decrease  indefinitely  for  a  pulsed  signal  of 
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bandwidth,  W ,  but  would  reach  a  minimum  width,  5Fminl  as  determined  by  the  transmitted  bandwidth 
to  radar  frequency  ratio: 

8  f mint/ Parang  =  SFBW/  FRadar  (16) 

where  FBraKK  is  the  Doppler  frequency  of  the  Bragg  line;  hFBW  is  the  transmit  bandwidth  used;  and 
FRujur  is  the  radar  operating  frequency.  For  a  100-/zs  pulse  (10-kHz  bandwidth),  and  a  10-MHz  radar 
frequency,  the  ratio  of  the  above  equation  is  0.001.  For  this  case,  the  minimum  line  width  one  can 
achieve  by  increasing  integration  time  is  0.001  times  the  Bragg  line  Doppler  shift,  or  0.00314  Hz,  where 
fB  (Hz)  =  0. 102^/,/radar  (MHz)  relates  the  radar  frequency  and  the  Bragg  line  Doppler  shift.  Narrow¬ 
ing  the  bandwidth  to  100  kHz  using  a  1 0-^t.s  pulse,  would  result  in  just  a  0.0314-Hz  minimum  achiev¬ 
able  Bragg  line  width.  Further  increases  in  the  integration  time  beyond  this  minimum  simply  creates 
more  Doppler  filters,  narrower  in  bandwidth,  the  sum  of  which  maintain  the  shape  of  this  minimum 
width  Bragg  line  signal.  This  is  the  minimum  one  should  achieve  if  the  scattering  surface  were  per¬ 
fectly  coherent  over  the  entire  radar  cell  (scatter  from  a  land  surface,  for  example)  and  if  the  iono¬ 
sphere  were  a  perfectly  coherent  reflector.  However,  earlier  width  limitations  are  reached,  which  are 
broader  than  Fmjnl. 

One  earlier  limitation  to  minimum  achievable  line  width  is  determined  by  patch  size  effects  on  the 
ocean  surface.  That  is,  the  scattering  surface  is  not  coherent  over  the  radar  pulse  length.  This  is  pri¬ 
marily  due  to  the  fact  that  the  ocean  wave  spectrum  is  not  a  linear  one,  but  does  exchange  energy 
between  components.  To  first  order,  the  Fourier  approximation  holds,  since  the  time  of  transfer  of 
energy  is  long  compared  to  most  measurements  of  a  given  component.  As  one  increases  the  integra¬ 
tion  time  on  processing,  one  reaches  times  of  the  order  of  energy  exchange  processes  within  the  wave 
spectrum.  The  waves  are  behaving  as  wave  packets,  equivalently,  with  spatial  extent  smaller  than  the 
radar  pulse.  Now  the  minimum  linewidth  one  can  achieve  for  a  given  Bragg  line  is  no  longer  related  to 
the  reciprocal  of  the  pulse  length,  but  to  the  reciprocal  of  the  size  of  the  wave  packets  on  the  ocean 
surface.  Although  this  effect  has  not  been  investigated  in  detail,  typical  values  of  100  to  200  s  are  the 
limitations  imposed  by  wave  packet  size  effects  to  the  integration  times  one  can  use  before  the  Bragg 
line  ceases  to  narrow. 

The  most  severe  linewidth  limitation  is  provided  by  the  ionosphere.  Propagation  through  the 
magneto-ionic  medium  and  refraction  back  to  earth  results  in  phase  changes  which  may  not  be  uniform 
across  the  phase  front  of  the  radar  wave.  If  one  thinks  of  the  refraction  as  a  reflection  from  a  virtual 
surface  at  some  virtual  height,  which  produces  the  equivalent  final  wave  vector,  one  can  again  think  in 
terms  of  coherent  areas  on  this  virtual  surface.  These  areas  are  due  to  ionospheric  turbulence,  and 
cause  time  rates  of  change  of  phase  paths  of  refracted  radar  signals.  For  F-layer  propagation  the  size  of 
these  patches  is  typically  less  than  those  on  the  ocean  surface,  and  the  maximum  useful  integration 
times  available  are  typically  25  s,  although  more  stable  conditions  are  often  encountered.  Propagation 
via  the  F-layer  is  generally  much  better,  and  quite  often  the  sea  surface  is  the  limiting  factor  in  deter¬ 
mining  minimum  achievable  line  widths.  Integration  times  of  up  to  100  s  have  been  used  successively 
via  stable  F-layer  propagation. 

The  net  effect  of  this  section  on  our  point  of  interest  is  that  one  must  consider  the  total  energy  within  the 
broadened  Doppler  Bragg  line  as  the  measurable  quantity  for  the  determination  of  both  p  and  C-  In  either 
the  case  of  the  sea  surface  or  the  ionosphere,  where  patch  size  effects  dominate,  scatter  from  the  vari¬ 
ous  patches  have  some  distribution  which  is  Gaussian  regarding  the  sum  phase  of  scatter  from  all  of  the 
resonant  patches.  In  processing  radar  data  for  these  cases  one  should  attempt  to  match  the  integration 
time  used  to  the  bandwidth  limiting  the  linewidth  so  that  the  peaks  of  the  Bragg  lines  are  a  good 
representation  for  the  energy  within  the  Bragg  line,  i.e.,  so  that  all  of  the  appropriate  energy  in  a  Bragg 
line  is  contained  in  a  single  Doppler  filter.  Then  one  can  simply  scale  peaks  in  measuring  {.  and  no 
additional  corrections  for  integration  times  or  linewidth  effects  are  necessary.  If  longer  integration 
times  are  used  than  this  "matched”  one,  one  must  sum  up  the  energy  in  all  the  Doppler  filters  which 
contain  significant  contributions  to  the  Bragg  line.  Shorter  integration  times  than  the  matched  one  may 
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be  used  with  Kttle  effect,  as  long  as  second-order  scatter,  sufficiently  far  from  zero  Doppler  so  that  Eq. 
(14)  is  poorly  satisfied,  does  not  fall  within  the  zero  Doppler  filter,  affecting  the  measure  of  (. 

An  example  of  the  effect  the  coherence  of  the  ionosphere  has  on  minimum  achievable  linewidth 
is  shown  in  Fig.  5.  A  total  time  period  of  409.6  s  of  radar  data  was  processed  with  several  integration 
times,  then  incoherently  averaged  so  that  the  same  data  period  is  covered  in  each  case.  For  the  transi¬ 
tion  from  25.6  to  51.2  s  occurs,  the  expected  line  narrowing  occurs,  and  a  3  dB  increase  of  the  Bragg 
peak  relative  to  the  continuum  second-order  occurs.  This  second-order  scatter,  which  is  broadband,  will 
fill  all  Doppler  bins  between  the  Bragg  lines  irrespective  of  the  coherent  integration  time  used  here.  As 
the  coherent  time  is  increased  from  51.2  to  102.4  s,  the  effects  of  the  ionosphere  begin  to  appear,  the 
line  does  not  narrow  by  a  half  as  expected,  nor  is  there  any  increase  in  the  ratio  of  Bragg  line  peak  to 
zero  Doppler  energy.  In  both  cases  now  the  scattered  energy  has  a  finite  inherent  line  width  imposed 
by  the  ionosphere.  Further  increases  in  integration  time  will  simply  fill  additional  Doppler  bins  with 
scattered  energy,  with  no  change  in  amplitude  or  improvement  in  signal  to  noise.  The  signal  has  in 
effect  become  noise-like  in  the  sense  that  it  is  spread  over  a  finite  bandwidth  and  cannot  be  further 
spectrally  isolated  by  increasing  integration  time.  The  final  change  from  102.4  to  204.8  s  again  causes 
no  change  in  the  ratio  of  Bragg  peak  to  zero  Doppler  continuum  level,  although  the  Bragg  lines  are 
decomposed  more  clearly  into  the  various  ionospheric  multipath  contributions. 


or0  DEPENDENCE  ON  SPREADING  FUNCTION  USED 

Determination  of  crfn  is  affected  in  two  major  ways  by  the  spreading  function  one  uses  to  define 
the  directional  sea  spectrum.  First,  the  value  of  p  is  determined  by  the  ratio  of  amplitudes  of  the  direc¬ 
tional  wave  spectrum  at  an  angle  parallel  to  the  direction  of  the  radar  look  divided  by  the  amplitude  n 
radians  to  that  direction,  F(O)/F(0  +  n).  Second,  the  amplitude  of  o-0  along  the  wind  direction  will 
vary  in  an  absolute  sense  for  different  angular  spreading  models,  and  these  are  known  to  differ  greatly 
for  early  periods  of  wave  growth,  when  the  wave  energy  is  highly  peaked  along  the  wind  direction, 
versus  later  periods  when  the  wave  energy  has  spread  in  angle.  (See,  for  example,  [14].)  In  addition  to 
this  fixed  frequency  behavior  with  time,  angular  spreading  is  observed  to  vary  with  wave  frequency  for 
a  given  time  in  the  wave  development  process  as  well.  As  a  result  of  this  type  of  behavior,  a  direc¬ 
tional  spectrum  highly  peaked  along  the  wind  direction,  compared  to  a  spectrum  generally  uniformly 
distributed  in  angle,  can  give  the  same  integrated  result  for  a  given  frequency  component  when  meas¬ 
ured  by  an  omnidirectional  instrument  such  as  a  Waverider  buoy.  However,  a  radar  measurement  will 
yield  a  greater  angular  spectral  power  value  for  a  look  angle  along  the  wind  direction  in  the  case  of  the 
highly  peaked  spectrum  than  for  the  spread  spectrum,  since  it  is  measuring  just  a  small  angular  portion 
of  the  ocean  wave  directional  spectrum  rather  than  the  integrated  value  over  all  angles  as  is  done  by  a 
buoy.  The  value  of  p  measured  will  differ  for  the  two  cases  as  well,  and  can  give  incorrect  results  in 
use  of  the  nomograph  if  the  sea  spectrum  does  not  spread  as  given  by  Eq.  (15).  However,  for  open 
ocean  conditions  and  far  from  the  spectrum  peak  in  frequency,  Eq.  (15)  can  be  considered  as  valid. 

Concerning  the  second  point,  Barrick’s  original  paper  defining  <r°  as  -29  dB  was  developed  using 
an  omnidirectional  wave  spectrum  with  equal  amplitudes  in  all  directions  in  the  forward  half  plane,  with 
no  energy  in  the  back  half  plane  opposite  the  wind,  and  with  a  constant  a  of  0.01  in  Eq.  (4).  For  this 
case,  p  could  not  be  defined  because  a  weaker  Bragg  line  does  not  exist.  If  one  sets  a  back  plane 
energy  level  some  constant  value,  some  number  of  dB  down  from  the  forward  plane  energy,  the  value 
of  p  remains  constant  with  look  angle,  which  observationally  is  not  true.  If  one  allows  some  direc¬ 
tionality  to  the  spreading  function,  then  the  amplitude  of  the  radar  cross  section  must  increase  along 
the  wind  direction  compared  to  the  -29  dB  predicted  by  Barrick,  such  that  the  integrated  wave  energy 
over  all  directions  is  the  same  for  any  directional  spreading  function. 
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Fig  5  —  The  same  period  of  radar  data,  409  6  s.  was  processed  using  four  dilTerem  in- 
legraiion  limes,  with  specira  incoherently  added  io  cover  an  identical  time  period  for 
skywavc  propagation.  With  the  change  from  25.6  to  51.2  s  < 5a > ,  narrowing  of  the  Bragg 
lines  and  the  expected  3  dB  increase  in  i  is  observed,  as  described  in  the  text.  With  a 
further  increase  to  102  4  <5b)  and  204.8  s  (5c).  ionospheric  multipath  prevents  further 
Bragg  line  narrowing,  and  no  change  in  y.  as  described  in  the  section  on  effects  of  in¬ 
tegration  time,  etc 


13 


D.  B.  TRIZNA 


3 


Figure  6  illustrates  this  fact  by  a  plot  of  o-0  as  a  function  of  angle  to  the  wind  for  a  selection  of 
models  of  directional  spreading  functions,  using  Phillips’  value  of  0.0081  for  a  in  place  of  Barrick’s 
0.01.  The  corresponding  constant  value  of  -29.9  is  the  smallest  value  of  the  set  along  the  wind  direc¬ 
tion,  while  the  largest  value  in  this  direction  is  given  by  the  most  highly  peaked  spreading  function.  In 
all  cases,  the  integrated  energy  under  each  curve  over  all  angles  is  the  same,  as  is  necessary  from  the 
discussion  in  conjunction  with  Eq.  (3). 


Fig.  6  —  Curves  of  <r°  for  "saturated"  sea  spectra  are  shown  as  a  function  of  radar-wind  angle  for 
different  models  of  the  directional  spreading  function.  F(0),  of  Eqs.  (3)  and  (4)  in  the  text.  The 
straight  line  is  Barrick's  result,  slightly  modified  for  at  Phillips  constant  of  0.0081  rather  than  0.005 
used  by  Barrick  [31.  It  is  seen  that  as  energy  is  localized  about  the  forward  direction  relative  to  the 
omnidirectional  model,  the  cross  section  along  the  wind  direction  must  rise  compared  to  the  con¬ 
stant  value  of  -29.9  dB.  This  normalization  comes  about  because  the  integrated  area  over  all  an¬ 
gles  under  the  curve  must  give  the  same  result,  corresponding  to  the  energy  measured  by  an  omni¬ 
directional  wave  measuring  device  which  does  not  distinguish  angular  information,  and  for  which 
the  Phillips,  Pierson-Moskowitz,  and  other  wave  models  are  defined. 


The  spreading  function  used  will  affect  the  determination  of  <r°  in  one  additional  way,  via  its  use 
in  calculating  <r&)  and  derivation  of  the  nomograph.  The  double  scatter  mechanism  used  to  calculate 
the  second-order  contribution  at  zero  Doppler  frequency  depends  upon  this  spreading,  but  the  magni¬ 
tude  of  this  effect  has  not  been  investigated.  An  example  of  the  variation  in  the  angular  spreading  of 
two  models  is  shown  in  Fig.  7,  where  we  have  plotted  the  spreading  function: 

ir~ 1  F(9)  -  /3  +  cos2#,  for  |#|  <  rr/2  (17a) 

7T_I  F(9)  -  /32/F(9  +  tr),  for  |#|  >  rr/2  (17b) 

as  well  as  an  example  of  another  popular  model  given  by: 

ir~'  F(9)  -  0  +  cos1  (0/2).  (18) 

Constants  0  were  chosen  for  both  models  so  that  the  front  to  back  ratio  are  both  24  dB.  The  exponent, 
s,  of  the  second  model  was  then  chosen  such  that  both  models  had  the  same  value  at  90°  to  the  wind 
direction  as  well.  It  is  noted  that  there  are  no  more  than  a  few  dB  variation  over  all  angles,  and  that 


14 


NRL.  REPORT  8579 


Fig.  7  —  Plots  of  the  spreading  function,  F(«),  of  Eqs.  (17)  and  (18)  of  the  text  are 
shown  plotted  as  a  function  of  angle  to  wind.  The  solid  curve  is  the  model  used  in 
Trizna,  et  al.  (9)  whereas  the  dashed  curve  is  of  a  type  determined  from  pitch  and 
roll  buoy  analyses.  Constants  have  been  chosen  for  this  model  so  that  both  models 
have  the  same  upwind-downwind  ratio,  24  dB,  as  well  as  equal  amplitudes  at  90°. 
Both  are  empirical  and  neither  has  been  derived  from  first  principles.  The  point  to 
note  is  that  both  models  can  be  mapped  into  the  other  with  just  a  few  dB  of 
difference  over  all  angles.  The  dashed  curve  used  to  generate  radar  spectra  and  the 
nomograph  would  probably  produce  results  similar  to  the  solid  curve. 


both  models  are  virtually  identical  out  to  60°  from  the  wind  direction.  It  is  not  expected  that  the 
nomographs  generated  using  these  two  models  would  differ  significantly  from  one  another,  although 
use  of  a  much  more  peaked  spectrum  would  produce  different  results.  Additional  work  in  this  area, 
using  frequency  dependent  spreading  functions,  is  ongoing.  Indeed,  accurate  representation  of  the 
spreading  of  wave  energy  as  a  function  of  wave  frequency  is  one  of  the  biggest  uncertainties  faced  in 
modeling  second-order  radar  Doppler  spectra. 

EXAMPLES  OF  RADAR  DATA  COMPARISONS 
WITH  WAVERIDER  BUOY  RECORDS 

The  technique  of  determining  <r°  from  radar  data  using  the  nomograph  was  applied  to  surface 
wave  data  collected  at  San  Clemente  Island  using  a  radar  with  a  15°  beamwidth  that  is  constant  with 
radar  frequency,  made  so  by  selectively  omitting  elements  with  increasing  radar  frequency.  The  radar 
was  run  at  ten  frequencies  and  two  bearings  separated  by  30°,  switching  frequencies  and  then  bearings 
on  a  pulse-to-pulse  basis,  so  that  in  effect  the  data  were  collected  simultaneously.  The  radar  spectra 
were  processed  with  a  200-s  integration  time,  that  for  which  the  nomograph  was  developed  so  that  no 
correction  is  required  for  the  integration  time.  The  ten  frequencies  used  run  from  slightly  above  2 
MHz  to  24  MHz.  Data  were  collected  for  three  range  bins,  sampled  time  delays  of  50,  100,  and  150 
/xs,  and  all  were  used  for  this  analysis,  p  was  measured  by  taking  the  dB  difference  of  the  peaks  of  the 
Bragg  lines  and  not  measuring  the  power  in  adjacent  filters.  This  produces  no  errors  as  long  as  there  is 
no  broadening  of  the  Bragg  lines,  which  there  did  not  appear  to  be  for  the  integration  time  used.  £ 
likewise  was  measured  by  scaling  from  the  strongest  Bragg  line  peak  to  the  level  near  zero  Doppler 
which  represented  an  average  power  level  in  this  region.  The  nomograph  was  then  used  to  calculate  A, 
the  number  of  dB  the  omnidirectional  sea  spectrum  is  down  from  the  traditional  Phillips  O-5  asymp¬ 
tote,  which  crosses  -33  dB  m2/Hz  at  a  one-hertz  frequency.  It  is  A  which  is  then  plotted  as  a  function 
of  radar  frequency  used  for  comparison  with  buoy  records.  The  Phillips  asymptote  for  a  -  0.0081  is 
drawn  on  the  buoy  records  of  Figs.  8-10,  as  well  as  lines  5  and  10  dB  down,  so  that  wave  power  spectral 
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RADAR  FREQUENCY  (MHz) 


Fig.  8  —  Measured  values  of  X  determined  from  radar  data,  using  the  nomo¬ 
graph  of  Fig.  4  are  plotted  versus  radar  frequency.  These  are  equivalently,  10 
log(awf<5/0.008l),  so  that  each  of  the  points  is  an  estimate  of  the  number  of 
decibels  the  ocean  wave  spectrum  is  down  from  the  traditional  Phillips  asymp¬ 
tote  at  the  Bragg  resonant  ocean  wave  frequency  (given  roughly  by 
fB  (Hz)  =  0. 1 02>//‘radar  (MHz).  Circles  denote  the  right  of  two  beams,  while 
crosses  denote  the  left  beam.  Radar  spectra  which  were  noisy  at  zero  Doppler 
were  not  used.  Below  these  data  is  an  ocean  wave  spectrum  taken  simultane¬ 
ously  within  the  scattering  area  using  a  Waverider  buoy.  Dotted  lines  are 
drawn  through  the  Bragg  resonant  wave  frequencies  for  each  radar  frequency 
running  left  to  right.  Wave  amplitude  differences  from  the  top  near-diagonal 
line  representing  the  Phillips  asymptote  («  =  0  0081)  are  to  be  compared  with 
radar-determined  X  for  each  frequency  from  left  to  right  in  corresponding  ord¬ 
er.  The  lower  two  diagonal  straight  lines  are  down  5  and  10  dB  from  the  Phil¬ 
lips  asymptote,  respectively. 
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Fig.  9  —  Data  for  a  second  day  for  which 
near  simultaneous  radar  and  buoy  data  are 
available  for  comparison,  for  a  case  closer  to 
saturation  than  that  in  Fig.  8. 
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Fig.  10  —  Data  sim to  8  and  9,  but  for  a  case  even  closer 
.o  tht  Philips  saturation  reference. 


value  differences  from  the  Phillips  asymptote  can  be  compared  directly  with  the  value  of  A  derived 
from  the  radar  data  for  each  Bragg  vave  frequency. 

Figure  8  shows  the  results  of  such  a  comparison  for  nine  frequencies,  three  of  them  grouped 
within  20  kHz  of  one  another  near  13.4  MHz.  Right  and  left  beam  data  are  plotted  separately,  and  are 
a  good  indication  of  the  self-consistency  of  the  analysis.  The  data  for  the  24-MHz  frequency  were  omit¬ 
ted  because  the  noise  level  was  higher  than  the  zero  Doppler  continuum  and  no  value  of  A  could  be 
determined.  This  was  also  the  case  for  the  second  and  third  range  bin  for  the  20-MHz  data,  due  to  sur¬ 
face  wave  propagation  losses  lowering  the  signal  relative  to  the  noise.  Dotted  lines  are  drawn  indicating 
the  Bragg  ocean  wave  frequencies  sensed  by  each  of  the  radar  frequencies  in  the  buoy  spectrum  below. 
The  buoy  was  in  the  center  of  one  of  the  radar  cells,  and  the  spectrum  was  taken  at  the  time  indicated. 
The  value  of  A  determined  from  each  of  the  two  beam  positions  is  indicated  by  different  symbols.  The 
general  overlap  of  the  two  types  of  symbols  indicates  that  the  spreading  function  model  used  in  deter¬ 
mining  the  nomograph  is  not  unreasonable,  since  both  beam  positions  give  similar  and  independent 
values  of  A.  The  average  value  of  A  for  each  radar  frequency  is  to  be  compared  with  the  number  of  dB 
the  corresponding  ocean  wave  power  level  is  down  from  the  Phillips  asymptote,  based  upon  a  value  of 
a  of  0.0081  in  the  theory.  Generally,  the  agreement  is  quite  good,  particularly  at  the  higher  frequen¬ 
cies  where  the  spreading  function  model  can  be  expected  to  hold  rather  well.  The  worst  disagreement 
is  for  the  lowest  radar  and  wave  frequency  and  is  probably  due  to  the  spreading  being  much  narrower 
than  that  assumed  in  the  model.  Note  the  complex  behavior  of  this  wave  spectrum  with  ocean  wave 
frequency,  a  very  good  example  of  a  case  which  could  not  be  fit  by  a  model  such  as  the  JONSWOP. 
Note  that  any  inversion  techniques  for  determining  wave  parameters  from  radar  parameters  which 
depend  upon  a  given  model  in  implementing  the  technique  would  have  problems  with  this  spectrum. 

A  second  example  is  shown  in  Fig.  9,  which  is  much  closer  to  the  Phillips  asymptote  than  the  pre¬ 
vious  case.  Again  data  were  not  available  for  the  24  MHz  frequency  and  some  data  points  are  omitted 
for  other  radar  frequencies  also,  because  of  noise  overwhelming  the  zero  Doppler  region.  Again  the 
agreement  is  quite  good  with  the  worst  case  being  the  lowest  radar  frequency. 
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The  third  example  shown  in  Fig.  10  is  for  the  case  closest  to  saturation  of  the  three.  Again  data 
were  omitted  when  noise  precluded  a  f  measurement.  The  agreement  between  the  two  data  sets  is 
again  quite  good. 

A  fourth  example  is  shown  in  Fig.  11  for  a  case  of  cross  wind  data  for  one  of  the  look  angles,  but 
for  which  there  were  no  wave  records  collected.  The  point  of  this  data  is  to  show  that  the  technique 
appears  to  predict  consistent  results  for  the  highest  frequencies,  where  both  sets  of  data  for  the  two 
look  angles  appear  distributed  similarly.  As  one  drops  in  radar  frequency,  the  results  for  the  two  beams 
tend  to  diverge,  indicating  again  that  the  angular  spreading  of  the  directional  spectrum  is  probably  nar¬ 
rower  for  the  lower  frequencies  than  is  described  by  theory,  and  may  affect  the  cross  wind  case  in  a 
more  severe  manner  than  other  cases  for  the  lower  frequencies. 


RADAR  FREQUENCY  (MHz) 


Fig.  1 1  —  Data  similar  to  the  last  three  figures,  but  for  which  no  buoy  spectra 
were  available  These  data  are  of  interest,  however,  because  of  the  fact  that 
the  crosses  represent  cases  for  radar  wind  angles  near  ninety  degrees.  It  gives 
some  estimate  of  where  the  nomograph  accuracy  becomes  poor,  as  determined 
from  the  poorer  self-consistency  of  the  data,  as  described  in  the  text. 


«r°  ESTIMATION  FROM  A  MEASUREMENT 

The  value  of  the  radar  cross  section  of  the  sea  surface  can  now  be  estimated  using  the  inferred 
value  of  A  from  the  preceeding  analysis.  It  is  simply  the  value  of  cr°  from  Fig.  6  for  the  cosine  squared 
plot  at  zero  degrees,  -26.9  dB,  less  A,  less  the  angular  term  due  to  the  spreading  function  used  in  the 
theory,  the  solid  curve  of  Fig.  7.  In  the  last  term  we  also  account  for  the  second  Bragg  line,  which  is 
important  for  near  cross  wind  conditions.  (Strictly  speaking,  one  should  also  include  the  second-order 
continuum,  since  a-0  is  defined  as  the  sum  over  all  Doppler  frequencies  of  energy  scattered  from  the 
sea.  However,  the  second-order  scatter  is  far  less  than  either  of  the  two  Bragg  first-order  terms  and  will 
be  neglected  here.)  One  can  then  write  an  equation  for  <r°: 

<r°<0)  -  <r°(0)  +  A  +  G(fl). 

where  the  model  of  Eq.  (17)  is  used  for  360  degree  spreading: 
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G{0)~  10  log  lF(0)  +  F(0  +  it)]  (20a) 

-  10  log  Fi.it/ 2)  +  10  log  [F{0)/ F(n/D  +  F(0  +  it)/ Fi.ir/2))  (20b) 

-  -  12  dB  +  p/2*  10  log  (1  +  F(0  +  it)/F(0)]  (20c) 

-  -  12 dB  +  p/ 2*10  log  [1  +  IOtp/'0].  (20d) 

Note  that  the  second  part  of  the  last  term  will  only  be  important  for  cases  near  crosswind  when  both 
Bragg  lines  contribute  nearly  equally.  Otherwise,  this  term  is  small  and  o0  will  be  determined  by  the 
amplitude  of  the  strongest  Bragg  line.  Using  the  value  of -26.9  dB  for  the  upwind  value  of  <r°  gives 

<r°  (0)  =  -  38.9  +  A  +  (p/ 2)*10  log  [1  +  KT"'10].  (21) 

APPLICATION  TO  ESTIMATION  OF  TARGET  CROSS  SECTION 

As  a  final  example  of  application  of  the  technique,  we  estimate  a-0  for  a  case  of  skywave  propaga¬ 
tion  in  which  a  target  is  seen,  most  probably  a  ship,  and  estimate  its  cross  section.  Figure  12  shows  the 
radar  spectrum,  collected  with  a  17-s  integration  time  and  a  21.8-MHz  radar  frequency,  at  a  range  of 
about  1600  km.  The  two  Bragg  lines  are  separated  by  just  a  few  dB,  indicating  near  cross  wind  condi¬ 
tions,  and  £  is  some  16  dB  below  the  strongest  Bragg  line.  From  the  nomograph,  £  for  this  case  is  24 
dB  for  the  200-s  integration  time.  The  correction  due  to  the  17-s  integration  time  is  the  log  of  the 
ratio,  17/200,  about  10  dB,  yielding  a  theoretical  £  of  14  dB,  which  is  2  dB  larger  than  the  measured 
value.  Hence  a  value  of  A  of  —2  dB  is  obtained. 


Fig.  12  —  An  example  spectrum  is  shown  in  which  a  ship  target  return  is  separated  from 
the  clutter,  with  its  cross  section  estimated  from  the  determination  of  A  and  tr°.  using 
the  technique  developed  here.  Note  that  these  values  for  it0  and  the  cross  sections  do 
not  include  the  additional  12  dB  some  authors  add  for  skywave  propagation.  By  our  con¬ 
vention,  skywave  and  surface  wave  cross  sections  are  identical,  and  the  12  dB  is  added  as 
a  propagation  gain  to  account  for  the  received  power  increase  of  12  dB  for  skywave  cases 


Using  the  results  derived  up  to  this  point,  we  now  estimate  tr0  for  the  example  of  Fig.  12,  where 
A  —  -2  dB  and  p  —  2  dB.  It  follows  from  Eq.  (21)  that,  for  A  *  -2  dB  and  p  —  2  dB: 

< r°(0 )--  38.9  -  2  +  1  *(2.12) -  38.8  dB  m2/m2.  (22) 

The  scattering  cell  for  the  pulse  length  and  beam  width  used  was  of  the  order  of  90  dB  above  a  square 
meter,  so  that  the  sea  surface  scattering  cross  section  was  51.2  dB  m2.  The  value  of  the  cross  section  of 
the  target  is  then  scaled  from  this  amplitude  of  the  greater  Bragg  peak,  which  is  22  dB  higher,  and  is 
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therefore  29.2  dB  above  a  square  meter,  or  832  m2.  Note  that  some  authors  add  12  dB  to  the  cross  section 
for  sky  wave  propagation  to  account  for  the  antenna  image  and  the  target  image.  We  choose  to  include  these 
as  propagation  effects  in  the  radar  equation. 

One  is  cautioned  in  applying  this  technique  to  target  cross  section  estimation,  to  be  very  careful  to 
choose  the  largest  value  for  the  target  signal  out  of  a  series  of  Doppler  spectra,  because  polarization 
fading  on  targets  is  far  more  severe  for  a  small  target  than  the  Doppler  spectrum  of  the  sea  surface. 
Ideally  one  should  average  the  Doppler  sea  echo  spectrum  for  a  minimum  of  10  min,  and  take  the  larg¬ 
est  target  cross  section  observed  in  that  time  period  for  comparison. 

SUMMARY 

We  have  developed  a  technique  that  allows  one  to  estimate  the  cross  section  of  the  sea  surface  at 
HF  by  a  very  simple  measurement  of  two  characteristics  of  the  radar  Doppler  spectrum:  the  approach- 
recede  Bragg  line  ratio,  p,  and  the  spectral  energy  ratio  of  the  strongest  Bragg  line  to  the  average  con¬ 
tinuum  value  near  zero  Doppler  frequency,  f .  Using  a  nomograph,  based  upon  a  theoretical  model  for 
the  Doppler  spectrum,  one  can  estimate  the  scattering  cross  section  per  unit  area,  <r°,  for  the  observed 
patch  of  sea  surface  being  illuminated.  The  technique  therefore  does  not  require  any  system  calibration 
or  knowledge  of  antenna  gain  patterns,  ionospheric  D-region  losses,  and  the  like,  and  so  is  self  calibrat¬ 
ing.  The  technique  was  applied  to  radar  data  collected  in  the  surface  wave  mode  and  the  results  were 
compared  to  independent  measurements  of  the  omnidirectional  sea  spectrum  made  by  a  Waverider 
buoy.  Since  the  data  agree  quite  well  with  one  another  for  the  several  look  angles  and  ranges  that  were 
used  for  the  measurements,  the  implication  is  that  the  model  for  the  spreading  function  used  to 
develop  the  nomograph  is  generally  a  valid  one  for  most  radar  frequencies.  The  data  showed  a  ten¬ 
dency  to  disagree  with  the  buoy  spectra  for  the  lowest  radar  frequencies,  near  2  and  4  MHz,  and  this  is 
thought  to  be  due  to  the  prevailing  directional  spreading  function  not  being  well  represented  by  the 
model  that  was  used  for  these  ocean  wave  frequencies.  The  wave  spectrum  is  known  to  be  very  narrow  for 
frequencies  near  the  ocean  wave  spectrum  peak,  as  well  as  for  sea  conditions  which  are  newly  developing,  and 
our  technique  should  be  applied  with  caution  for  these  cases.  To  date,  the  uncertainty  in  knowledge  of  the 
variability  of  the  spreading  function  with  frequency  and  changing  sea  conditions  appears  to  be  the  big¬ 
gest  uncertainty  in  understanding  the  Doppler  spectrum  at  HF.  This  work  is  meant  to  be  an  introduc¬ 
tion  to  using  the  radar  Doppler  spectrum  as  a  cross  section  calibration,  and  it  is  felt  to  be  sufficiently 
accurate  to  be  applied  to  most  conditions  appropriate  to  ionospheric  propagation,  in  particular  the  high 
frequency  portion  of  the  HF  band.  Work  will  continue  in  updating  and  refining  techniques  introduced 
here. 
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